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A STABLE MARRIAGE OF POISSON AND LEBESGUE 

By Christopher Hoffman, 1 Alexander E. Holroyd 2 
and Yuval Peres 3 

University of Washington, University of British Columbia and UC Berkeley 

Let 3 be a discrete set in R d . Call the elements of S centers. The 
well-known Voronoi tessellation partitions R d into polyhedral regions 
(of varying sizes) by allocating each site of R d to the closest center. 
Here we study "fair" allocations of R d to 3 in which the regions 
allocated to different centers have equal volumes. 

We prove that if 3 is obtained from a translation-invariant point 
process, then there is a unique fair allocation which is stable in the 
sense of the Gale-Shapley marriage problem. (I.e., sites and centers 
both prefer to be allocated as close as possible, and an allocation is 
said to be unstable if some site and center both prefer each other over 
their current allocations.) 

We show that the region allocated to each center £ is a union 
of finitely many bounded connected sets. However, in the case of a 
Poisson process, an infinite volume of sites are allocated to centers 
further away than £. We prove power law lower bounds on the allo- 
cation distance of a typical site. It is an open problem to prove any 
upper bound in d> 1. 

1. Introduction. Let d> 1. We call the elements of R d sites. We write | • | 
for the Euclidean norm and C for Lebesgue measure or volume on M. d . Let 
3 C M. d be a discrete set. We call the elements of 5 centers. Let a € [0, oo] 
be a parameter, called the appetite. An allocation (of M. d to 3 with appetite 
a) is a measurable function 

^:R d ^SU{oo,A} 



Received June 2005; revised October 2005. 

Supported in part by MSRI, NSF Grant DMS-00-99814 and the University of Wash- 
ington Research Royalty Fund. 

2 Supported in part by an NSERC (Canada) research grant and by MSRI, CPAM and 
PIMS. 

Supported in part by NSF Grants DMS-01-04073 and DMS-02-44479 and by MSRI 
and CPAM. 

AMS 2000 subject classification. 60D05. 

Key words and phrases. Stable marriage, point process, phase transition. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Probability. 
2006, Vol. 34, No. 4, 1241-1272. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



C. HOFFMAN, A. E. HOLROYD AND Y. PERES 



such that £[^~ 1 ( A )] = °) and 

for all £ € 3. We call V ;_1 (£) the territory of the center £. We say that £ is 
sated if £[V'~ 1 (£)] = a an d unsated otherwise. We say that a site x is claimed 
if ^(cc) £ E and unclaimed if ?/>(x) = oo. Note that oo and A will have very 
different meanings: ip{x) = oo means that x is unable to find any center 
willing to accept it, whereas ij)(x) = A means that tp(x) is "undefined" — we 
allow an £-null set of such sites purely as a technical convenience. 

The following definition is an adaptation of that introduced by Gale and 
Shapley [6]. The idea is that sites and centers both prefer to be allocated as 
close as possible, and an allocation is said to be unstable if some site and 
center both prefer each other over their current allocations. 

Definition of stability. Let £ be a center and let x be a site with ip(x) £ 
{£, A}. We say that x desires £ if 

\x — £| < \x — ip(x)\ or x is unclaimed. 

We say that £ covets x if 

\x — £| < \x' — £| for some x' 6 2/> _1 (£), or £ is unsated. 

We say that a site-center pair (a?,£) is unstable for the allocation ip if x 
desires £ and £ covets x. An allocation is stable if there are no unstable pairs. 
Note that no stable allocation may have both unclaimed sites and unsated 
centers. 

Figures 1 and 2 illustrate stable allocations for several values of a on the 
2-torus. 

We consider three main questions about stable allocations. First, given a 
countable set in the plane, does there exist a stable allocation? Second, if 
there is a stable allocation, is it unique? Finally, what do the territories look 
like for a typical set of centers, qualitatively and quantitatively? 

Theorem 1 (Existence). For any discrete set of centers S C M. d , and 
any appetite a€ [0, oo], there exists a stable allocation. 

Note that some condition on S is needed to guarantee existence of a stable 
allocation. For example, if d = 1 and H = {1, 1/2, 1/3, . . .}, then it is easy to 
see that some center must form an unstable pair with the site —1. 

We can establish uniqueness of the stable allocation if H is a sufficiently 
well-behaved set in two dimensions, or if 3 arises from a point process. 
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Fig. 1. A stable allocation in the critical case a = 1. (If you are looking at greyscale 
images, color versions are available at www.math.ubc.ca/~holroyd or math. PR/0505668.) 

Theorem 2 (Uniqueness for recurrent configurations). Given discrete 
H C let Gb be the graph with vertex set H and with an edge between every 
pair of centers at distance at most b from each other. Suppose S is such that, 
for all b > 0, the graph Gb is recurrent for simple symmetric random walk. 
Then for any a € (0, oo), there exists an C-a.e. unique stable allocation. 

Now let II be a translation-invariant simple point process on M. d , with 
intensity (i.e., mean points per unit volume) A 6 (0,oo) and law P. The 
support of IT is the random set [II] = {z G M d :LI({z}) = 1}. We consider 
stable allocations of the random set of centers S = [II] . 
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Fig. 2. Stable allocations with appetites a = 0.2, 0.8, 1, 1.2, 2, oo. The centers are cho- 
sen uniformly at random in a 2-torus, with one center per unit area. Each territory is 
represented by concentric annuli in two colors. 
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Theorem 3 (Almost sure uniqueness). For a translation-invariant point 
process II of finite intensity on W 1 and for any a £ (0, oo), there exists P-a.s. 
an C-a.e. unique stable allocation to S= [II]. Furthermore, can be 
chosen to be a measurable, isometry-equivariant factor of II. 

Note that uniqueness can fail for a general set of centers, at least in 
the setting of sufficiently general spaces. For example, let d=l, and E = 
{X}£=i l/i'n = 1,3,5, . . .}, and consider instead of C the measure with a unit 
point mass at each of Ya=i l/i'-n = 2,4,6, ... . With a = 1, there are multiple 
stable allocations: one assigns each mass of sites to the center immediately on 
its left; another assigns them to the right and leaves the center 1 unsated. We 
do not know whether uniqueness can fail for stable allocations of Lebesgue 
measure to a set of centers in W 1 . 

Henceforth, we let \& = be as in Theorem 3. 

Theorem 4 (Phases). Consider a stable allocation with appetite a to a 
point process of intensity A € (0, oo) which is ergodic under translations. 

(i) If Xa < 1 (subcritical), then a.s. all centers are sated but there is an 
infinite volume of unclaimed sites. 

(ii) If Xa = 1 (critical), then a.s. all centers are sated and L-a.a. sites 
are claimed. 

(hi) If Xa > 1 (supercritical), then a.s. not all centers are sated but L-a.a. 
sites are claimed. 

Note that by scaling M. d there is no loss of generality in assuming A = 1 
or a = 1. [The assumption of ergodicity can of course be weakened to cover 
any mixture of ergodic point processes, all of whose intensities satisfy the 
appropriate condition in (i)-(iii).] 

Note that the assertions of Theorem 4 may fail when generalized in the 
most obvious way to a general set of centers S. Indeed, let d = 1 and 
a = 1 and let E = Z \ {0}. The set E has "intensity" 1 in the sense that 
#{E n [— r, r]} ~ C[—r,r\ as r — > oo, so we might expect the configuration 
to be "critical." However, the stable allocation has a.e. site in (—1/2,1/2) 
unclaimed. (This may be deduced from a monotonicity result, Proposition 
21, by comparison with the case S = Z.) 

We now consider the geometry of territories. 

Theorem 5 (Geometry of territories). Consider a stable allocation to a 
translation-invariant point process of finite intensity. 

(i) All territories are bounded a.s. 

(ii) Any bounded subset of M, d intersects only finitely many territories 

a.s. 
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(iii) rnay be chosen to be a stable isometry-equivariant factor of H in 
such a way that a.s. each territory is a union of finitely many open connected 
sets. 

We now restrict to the case when II is a homogeneous Poisson process 
with intensity 1. The following result contrasts with Theorem 5. 

Theorem 6 (Infinite demand). Let II be a Poisson process and consider 
the critical case A = a = 1. 

(i) A.s. every center is desired by an infinite volume of sites. 

(ii) A.s. every site is coveted by infinitely many centers. 

[Note that (i) above holds trivially in the subcritical case, while (ii) holds 
trivially in the supercritical case, by Theorem 4.] 

We now consider quantitative properties of stable allocations. The quan- 
titative behavior turns out to be very different in the critical and noncritical 
phases. Broadly speaking, we find exponential tails in noncritical cases and 
power law tails in critical cases. A natural quantity to consider is 

X = |* n (0)| 

(i.e., the distance from the origin to its center), where we take |oo| = oo. Let 
E be the expectation operator associated with P. 

Theorem 7 (Critical lower bounds). Let II be a Poisson process with 
intensity A = 1. 

(i) For d=l,2 and a=l, we have EX d / 2 = oo. 

(ii) For d>3 and a = l, we have EX d = oo. 

Proofs of the following two theorems appear in a second article by the 
same authors [8]. 

Theorem II- 1 (Critical upper bound). Let II be a Poisson process with 
intensity A = 1. For d = 1 and a = 1, we have EX 1 / 18 < oo. 

We have not been able to prove any quantitative upper bound for X when 
d>2. 

Theorem II-2 (Noncritical upper bounds). Let IT be a Poisson process 
with intensity A = 1 . 

(i) For all d and a > 1, we have ~Eie cXd < oo; 

(ii) for all d and a < 1, we have ~Ei(e cX ;X < oo) < oo, 
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where c = c(d, a) > 0. 

It is straightforward to see that the bounds in Theorem II-2 are tight up 
to the value of c, since the distance from the origin to the closest center has 
tail decaying exponentially in volume. Some of the above bounds may be 
strengthened to bounds involving the radius of a typical territory; see [8] for 
more details. 

Remarks, (i) The concept of stability was introduced by Gale and 
Shapley [6] in the context of discrete "marriage problems," and has since 
been studied extensively; see [12] for more information. The model consid- 
ered here appears to be the first extension to a continuum setting. 

(ii) The cases a = 0, oo are relatively simple. When a = 0, every center's 
territory is either empty or consists only of the center itself. When a = oo, 
every center covets every site, so each tp(x) equals the closest center to x and 
we obtain the Voronoi tessellation (see [14] for more information). Thus, we 
restrict our attention to the cases a € (0, oo). 

(iii) A stable allocation has the following informal interpretation. De- 
note the ball B(x,r) = {y € R rf : \x — y \ < r}. Fix S and a, and imagine a 
sphere centered at each center, simultaneously growing linearly, so that, at 
time t, all spheres have radius t. Now let each center capture all sites en- 
countered by its sphere, provided the center is not yet sated, and provided 
no other center captured the site first. It is easy to convince oneself that if 
this informal model exits, it should give rise to a stable allocation, and the 
territory captured by £ up to time t should be V ; " 1 (C) nJB(£,t). However, 
it is by no means clear a priori that the model does exist, since the pro- 
posed evolution involves long-range dependence. (See [1] for an example of 
a seemingly natural model which turns out not to exist in certain settings.) 
One consequence of proving existence and uniqueness of stable allocations 
is that the informal picture described does indeed make sense. 

(iv) Our proof of existence (Theorem 1) will be a generalization of a 
proof in [10], which in turn is an adaptation of the algorithm introduced in 
[6]. In [10], stable allocations were used to construct certain shift-couplings 
of point processes. 

(v) It should be noted that stable marriage problems do not in general 
have unique solutions; see [6]. In this case the key to uniqueness (Theorems 
2 and 3) is that the preferences of sites and centers are "consistent," both 
being based on Euclidean distance. 

(vi) The condition on S in Theorem 2 is natural when d < 2. In par- 
ticular, by Corollary (3.52) in [16], it holds provided, for any b < oo, we 
have 

#{{£,£'} CEn.B(0,r): < b} = 0(r 2 log r) as r -> oo. 
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(vii) In the course of proving our main theorems, we will obtain several 
other results which are interesting in their own right, and which may be 
summarized as follows. Stable allocations are monotonic in both E and a in 
the sense that adding extra centers or increasing the appetite makes sites 
happier but centers more unhappy (Propositions 21 and 22). Although stable 
allocations are unique only almost everywhere, there is a unique "canoni- 
cal version" having the properties in Theorem 5 and other good properties 
(Theorem 24). In the second article [8] we prove that stable allocations are 
robust to far-away changes in S, in the sense that if a sequence of sets of 
centers converges weakly in the discrete topology, then the sequence of sta- 
ble allocations converges almost everywhere. The supercritical phase is rigid 
in the sense that if a is only just greater than one, a large proportion of 
centers are sated. 

(viii) Note that there is a partial symmetry between the roles of sites and 
centers in the definition of stability. In light of this, it is not surprising that 
many of our results (and the arguments used to prove them) come in "dual" 
pairs obtained by exchanging the roles played by sites and centers. 

2. Existence. In order to prove Theorem 1, we shall construct an allo- 
cation tfj using a variant of the Gale-Shapley stable marriage algorithm [6]. 
The algorithm has two different versions, called "site-optimal" and "center- 
optimal." (The reason for the names will become apparent in Lemma 15.) 
For our present purpose, one version will be sufficient, but later in the proof 
of Theorem 3 we shall have occasion to consider both. 

Site-optimal Gale-Shapley algorithm. Let W be the set of all sites which 
are equidistant from two or more centers. Since E is countable, W is £-null. 
For convenience, we take ip{x) = A for all x € W. 

We now construct ip by means of a sequence of stages. For each positive 
integer n, stage n consists of two parts as follows: 

(a) Each site x ^ W applies to the closest center to x which has not 
rejected x at any earlier stage. 

(b) For each center £, let A n (£) be the set of sites which applied to £ in 
stage n (a), and define the rejection radius 

r n (0 = inf {r : C(MO H r)) > a}, 

where the infimum of the empty set is taken to be oo. Then £ shortlists all 
sites in A n {^) n B(£, r n (£)), and rejects all sites in A n (^) \ B(£, r n (£)). 

We now describe tp. Consider a site x £ W. Since S is discrete, the fol- 
lowing is clear. Either x is rejected by every center (in order of increasing 
distance from x) or for some center £ and some stage n, x is shortlisted by £ 
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at all stages after n. In the former case we put ip(x) = oo (so x is unclaimed); 
in the latter case we put ip{x) = £. 

Proof of Theorem 1. Let be constructed according to the site- 
optimal Gale-Shapley algorithm above. We first check that ip is an alloca- 
tion. Let 5 n (£) be the set of sites shortlisted by a center £ at stage n. By the 
construction in (b) and the intermediate value theorem, we have C(S n (£)) < 
a. But by the definition of ip above, we have V ;_1 (£) = nmsu Prwoo Sn(€) = 
liminf n ^oo S n (£), so, by Fatou's lemma, we have < ot, as required. 

Note also that if a center £ ever rejects sites (at stage n say), then we must 
have £(S m (£)) = a for all stages m > n. Hence, an unsated center never 
rejected any sites. 

We must check that tjj is stable. Consider a site x and a center £. 

Suppose x desires £. If x is unclaimed, then it was rejected by all centers. 
On the other hand, if x is claimed and \x — £| < \x — tp(x)\, then the site x 
applied to ip(x) at some stage, therefore, it was rejected by the closer center 
£ at some earlier stage. In either case £ rejected x. 

Now suppose £ covets x. If \x — £| < \x' — £| for some x' £ ^ _1 (£), then 
£ never rejected x', so x' was never outside £'s rejection radius, and the same 
holds for the closer site x. On the other hand, if £ is unsated, then it never 
rejected any sites. In either case £ did not reject x. 

We deduce as required that (x,£) cannot be an unstable pair. □ 

3. Uniqueness for recurrent configurations. In this section we prove 
Theorem 2. We therefore consider a set 3 satisfying the assumption of the 
theorem. Without loss of generality (by scaling M. d ), we may also assume 
that a = 1. We work by contradiction. Suppose that tpi and V>2 are two sta- 
ble allocations to 3 which differ on an £-nonnull set. We will use them to 
define a certain flow on a directed graph with vertex set 3, and then show 
that this flow has properties which are impossible when 3 satisfies the given 
condition. 

For each center £, we start by defining several sets. Let the set of disputed 
sites for £ be the symmetric set difference 

Let 

where 

r(0 = sup{r :£[(^H0 U V>2 n J3(e, r)] < 1}- 
This is the set of disputed sites that £ prefers. Also define 

M(f ) = £>(£) n {x : \x - f | < max[|x - \x - ^2{x)\}}. 



10 C. HOFFMAN, A. E. HOLROYD AND Y. PERES 

This is the set of disputed sites that prefer £. 
Define for every pair of centers 

F(£,O = tf(0riA*(O- 

Let G be the directed graph with vertex set 3 and a directed edge from 
£ to £' whenever j£F(£,£') > 0. We label the directed edge (£,£') with the 
quantity £F(£,£'). We will regard these labels as capacities. We first show 
that G has no directed cycles. 

Lemma 8. For each center £, 

CF(M) = 0- 
Thus, G has no cycles of length one (self-loops). 

Proof. Suppose that £F(£,£) > 0. Then, without loss of generality, 
there exist a region A of positive volume and a center £' ^ £ such that: 

(i) Acj\r(0nM(0, 

(ii) ACiPY 1 ^) and 

(iii) AC^(e'). 

We claim that either 

(a) CWtHO \ B(£,r(0)) > CA (> 0), or 

(b) £(^(e))<l. 

To see this, suppose, contrary to (b), that Cip^iO = 1- The definition of 
r(£) then implies 

i = C([^(0 u ^HO] n r(0)) 

= ^(^(e) n s(£, r(0)) + ^r 1 ® \ ^ 2 " 1 (0] n r®)). 

Conditions (i), (ii) and (iii) above imply A C [V^CO \ ^2 HO] n B ((,, r (0), 
therefore, n B(£,r(£))] < 1 - £A Since Ctp^iO = 1. this implies 

(a) above. 

Now consider a site x£ A. We have that x G iV(£) C f?(£, r(£)). In case 
(a), for any y G t/VHO \ B (£>K£))> we have 

Hence, in either case (a) or (b), £ covets x under tp2- Also, since x G M(£), 
we have that 

\x-£\ < \x-?\, 

therefore, x desires £ under ip2- Thus, is an unstable pair for ip2, a 

contradiction. □ 
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Lemma 9 (Monotonicity of distances). For centers £1 7^ £2 and x\ G 
F{£,i->£,2), we have the following: 

(i) |£i - xi\ > \x\ - £ 2 |, and 

(ii) for any center £ and xq G F(£o, £1) \ iV(£i), 

I sco -£1 1 > |£i 

[Note that, by Lemma 8, the condition in (ii) above applies to a.e. site 
z eF(£ ,&).] " 

Proof. Since xi G M(£2), statement (i) follows from the definition of 
M(£ 2 ). For (ii), since x\ G iV(£i) and x G £>(£i) \ iV(£i), the definition of 
iV(£i) implies that |a?i — £i| < t(£i) < \x§ — £i|. □ 

Lemma 10. G has no directed cycles. 

Proof. By Lemma 8, there are no cycles of length 1. Suppose there is a 
cycle with vertices £o,£i,£ 2 , • ■ • ,£n = £o- Then for all i G {0, ... ,n — 1}, there 
exists Xi G F(£j, £j + i). By Lemma 8, we can further require that X{ £ ./V(£j+i). 
Then Lemma 9 yields 

|£o ~x \ > \x -£i| > |£i -rci| > ••■ > \x n - X - £ n | > |£ n - cc | 

Lemma 11 (Increasing capacity). For any center £, we /iaue 
^£F(£',£)< £ £F(£,£"). 

PROOF. By the definitions of F(£, £'), ]V(£) and M(£), we have that 

(1) [J F(£',£)CM(£) 

and 

(2) u ma=N(o- 

(Here U denotes disjoint union.) We claim now that 

£M(£)<£iV(£), 

which combined with (1) and (2) proves the lemma. To prove the claim, we 
consider two cases. First, if 

£iV(£)<i-£(Vr 1 (0nV2" 1 (6), 
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then r(£) = oo, so N(£) = D(g), and by Lemma 8, it follows that £M(f) = 0, 
establishing the claim. Second, suppose 

(3) £N(t) = l-£(^\0nip2 1 (0). 
By Lemma 8, we have 

2£(vr 1 (0 n v 2 _1 (0) +^(0 +>cm(o < ^vrHo+^Ho < 2, 

and combining this with (3) yields 

£M(e)<i-£(vr 1 (6n^ 2 - 1 (e)) J 

again establishing the claim. □ 

Supposing that ifti and ^2 do not agree £-a.e., there exists some center £ 
with £D(Q > 0, and therefore, CN(() > 0. Any site x € D(C) lies in M(C') 
for some center Therefore, G has some edge (C>C0* Now fix a directed 
edge (£, CO °f ^ an d let H be the directed subgraph consisting of this edge 
together with all vertices and edges of G that can be reached by directed 
paths from 

Lemma 12 (Bounded edges and summable paths). Let H be a graph as 
defined above. 

(i) The edges of H have uniformly bounded length. 

(ii) There exists K < oo such that, for any directed path (£o>£i)£2, ■ • ■) ™ 
H , we have 

oo 
i=0 

Proof. We start by proving (ii). Let (£o,£i>£2j . . .) be a directed path. 
Without loss of generality, we may assume that £o = C an d £i = C'- Define 

dj = essinf 16 — z\ 

and 

hi = ess sup |6 — z\. 

Clearly, bi > ai for each i, and we claim that also 

(4) ai>b i+ i, * = 0, 1, — 

To prove this, note that, by Lemma 9 (and Lemma 8), for £-a.e. Xj 6 
^(6>6+i) and every x i+ i G F(6+i,6+2), we have 

(5) |6 - xi| > |»i - 6+1 1 > I6+1 - Xi+i\; 
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then use the definitions of ai,bi + i. 

Let c = c(d) < oo be the volume of the unit ball. Using (4), we have 

OO 00 

]T 6+i) < CD(^ ) + Y,(cbf - caf) 

i=0 i=l 

oo 

<2 + c^(af_ 1 -af) 

i=l 

= 2 + ca d . 

But ao = essinf zg ^(£ )( */) [£ — z\, and this quantity is finite and independent 
of the choice of path, so we obtain (ii). 

Turning to (i), note that for a directed path as above, for i > 1 and £-a.e. 
Xi € F(6,6+i)> we have from (5) and (4) that 

|6 - 6+1 1 < |& - Zil + |xj - 6+1 1 

<2|6-Xi| 
< 26i < 2a . 

Since every edge of lies in some directed path starting with the edge 
(CCO, (i) follows. □ 

We will now deduce Theorem 2 from the following result on random walks. 
The proof of Proposition 13 uses standard technology, and is deferred until 
later in this section. 

Proposition 13 (Transience from capacities). Let H = (V,E) be a lo- 
cally finite, directed graph with no directed cycles, and with a distinguished 
vertex pEV. Suppose there exists a 11 capacity" function q: E — > [0, oo) which 
satisfies the following three conditions: 

(i) For every vertex v ^ p, we have 

Y q( u > v )< Y q(v,w). 

u:(u,v)£E w. (v,w)£E 

(") Y q( u 'P) = ° and Y i{pi w )>°- 

u:(u,p)£E w : (p,w)£E 

(hi) There exists K < oo such that, for every directed path (uq,u\, . . .) in 
H , starting at uq = p, we have 

oo 

Y,q(. u i, u i+l) < K - 
i=0 
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Then simple symmetric random walk on the undirected version of H starting 
at p is transient. 

Proof of Theorem 2. Suppose S has two stable allocations which 
differ on an £-nonnull set, and construct the graph H as described before 
Lemma 12 above. By the construction of H and Lemmas 10, 11 and 12(H), H 
satisfies the assumptions of Proposition 13 with (/(£,£') = CF(£,£'), hence, 
it is transient (for simple symmetric random walk on the undirected version 
of the graph). However, by Lemma 12(i), the assumption on E and the fact 
that a subgraph of a recurrent graph is recurrent (Corollary (3.48) of [16]), 
H is recurrent. This is a contradiction. □ 

Proof of Proposition 13. Without loss of generality, we may assume 
q is normalized so that J2w. (p,w)eE 9(/°> w) = 1. By a unit flow on H, we mean 
a function / : E — > [0, oo) satisfying the following: 

(i) for every vertex v ^ p, we have 

u:(u,v)aE w: (v,w)£E 

(ii) £ f(u, P ) = and £ f(p,w) = l. 

u:(u,p)£E w:(p,w)eE 

First we claim that there exists a unit flow / with / <q on E. This is 
proved using the max-flow/min-cut theorem (see, e.g., [15], Theorem 2.19). 
As a notational convenience we extend q to be on (V x V) \ E. For any 
finite set of vertices S 9 p, conditions (ii) and (Hi) of the proposition together 
with the normalization assumption yield 

Y Y <l( v i w ) >EE feO, ™) - q{w,v)] > 1, 

so the claim follows. Now from condition (iv) of the proposition, we have 
also 

oo 

(6) if 

i=0 

for every directed path (uo,u\, . . .) starting at p. 

If (uq,ux, . . .) is a directed path starting at p, then there is a unit flow 
which takes the value 1 on each edge of the path and elsewhere. Let T 
be the set of all unit flows of this form. By Proposition 2.20 of [15], for any 
acyclic unit flow /, there exists a unit flow /o satisfying /o < / on E, and a 
probability measure /ion T, such that, for all e G E, 

(7) fo(e)=J^(e)dp( 7 ). 
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Now, apply this to the unit flow / constructed earlier, which is acyclic 
since / < q. The summability condition on paths (6) implies that, for any 
7 G T, we have 



So /o is a unit flow of finite energy, and it follows (see, e.g., [4]) that H is 
transient as required. □ 

4. Almost sure uniqueness. In order to prove Theorem 3, we introduce a 
second version of the Gale-Shapley algorithm which constructs a potentially 
different allocation. Our approach to proving uniqueness will be to show that 
the two Gale-Shapley allocations are extremal (in an appropriate sense), and 
then to show that they in fact coincide. 

Center- optimal Gale-Shapley algorithm. Let W be the set of all sites 
which are equidistant from two or more centers, and take ip(x) = A for all 
x 6 W . For each positive integer n, stage n consists of two parts as follows: 

(a) For each center £, let i? n (£) be the set of sites which have rejected £ 
at some earlier stage, and define the application radius 



The center £ applies to every site x £ f?(£,a n (£)) \ R n (£,) \ W. 
(b) Each site x ^ W shortlists the closest center of those which applied 
to x in stage n (a) (if any), and rejects any others which applied. 

We now describe ip. Consider a site x ^ W. If no center ever applies to x, 
we put ip(x) = oo. Otherwise, if x shortlists some center at some stage, then 
x will shortlist this or a closer center at all later stages. Since E is discrete, it 
follows that, for some center £ and some stage n, x shortlists £ at all stages 
after n. In this case we put ij)(x) = £. 

Lemma 14. For any 5, the center- optimal Gale-Shapley algorithm yields 
a stable allocation. 

Proof. We mimic the proof of Theorem 1 in the site-optimal case. 
To see that ip is an allocation, note that the set A n (^) of sites which £ 
applies to at stage n satisfies C{A n ) < a, and V' _1 (C) = ummi n->oo^4.n(£) = 
limsup n ^ DO A n (^), and apply Fatou's lemma. 



£/ (e) 7 (e)<K 



Therefore, by (7) and Fubini's theorem, 




an® = mf{a : C{B{£, a) \ R n (C)) > a}. 
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To check stability, consider a pair (x,£). If a; desires £, then £ did not 
apply to x, while if £ covets x, then £ did apply to x. Hence, (x,£) is not 
unstable. □ 

For any allocation ip and any r 6 [0, oo], define the following quantities. 
For any site x ^ V -1 (A)> 

9x{ip,r) = l[\x - ip{x)\ < r]; 

and for any center £, 

As r varies, the above quantities measure how good the allocation is from 
the point of view of a particular site or center respectively. 

The following result states that the site-optimal Gale-Shapley allocation 
is the best for sites and the worst for centers, and vice versa for the center- 
optimal allocation. (The analogous fact for discrete stable marriage problems 
is well known: the algorithm in which boys propose gives the best possible 
stable system of marriages for boys and the worst possible for girls — see 
[6, 12].) 

Lemma 15 (Optimality). Fix a set of centers H, let ip be any stable allo- 
cation, and let ■0 GS ,'0 GSc be respectively the site-optimal and center- optimal 
Gale-Shapley allocations to the same H. For all r € [0, oo], we have the fol- 
lowing: 

(i) for C-a.e. site x, 

g x {i> GS ,r)>g x (^,r)>g x (^ GSc ,r); 

(ii) for every center 

7C (V- GSc , r ) > 7? (ip, r ) > ^ (^ GS , r) . 

The proof of Lemma 15 is deferred until the end of this section. 

Stable allocations to a point process. Now let II be a translation-invariant 
simple point process on M. d with intensity A € (0, oo), law P and expectation 
operator E. The assumption of finite intensity ensures that [II] is discrete 
a.s. A factor allocation is a measurable map ^, which assigns to P-a.e. point 
measure tt an allocation of [n], and which is translation-equivariant in 
the sense that, for any translation T of if *& n (x) = £, then ^Ttt(Tx) = T£ 
(here it is understood that ^t-k is defined if and only if ^ is). A factor al- 
location is isometry-equivariant if the above holds for all isometries T. Note 
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that the two Gale-Shapley algorithms clearly yield isometry-equivariant fac- 
tor allocations; we denote them \&^ s , *I'^ 1Sc , respectively. For a factor alloca- 
tion we shall consider the random allocation ^n> and when there is no 
risk of confusion, we denote it simply \&. Note that the joint law of the pair 
(IT,^n) is then translation- invariant. 

Let II* be the Palm version of the point process II, with law P* and 
expectation operator E*. (Recall that II* may be thought of as II conditioned 
to have a center at the origin, and that when II is a Poisson process, II* is a 
Poisson process with an added center at the origin; see [11], Chapter 11 for 
details.) By results of Thorisson [17] (see also [10] and [11]), II* and II may 
be coupled so that one is almost surely a translation of the other. Hence, if 
ty m is a factor allocation, then ^n* is defined P*-a.s. 

Lemma 16 (Site-center consistency). //II is a translation-invariant point 
process of intensity A £ (0, oo) and is a factor allocation of U, then 



The key to Lemma 16 is that the preferences of sites and centers expressed 
in the functions 5,7 are "consistent" — both being defined in terms of Eu- 
clidean distance, so that on average they coincide. The proof is deferred to 
the end of this section. 

Proof of Theorem 3. Note that YpGS^ $GSc are c j ear iy isometry-equi- 
variant factor allocations. Taking expectations in Lemma 15 and applying 
Lemma 16, we obtain 



Hence, we have equality throughout. Hence, by translation-invariance, we 
have 



wherever both quantities are defined. Hence, for all sites x we have P-a.s. 



E5o(^n,r) = AE* 7o (*n*,r 



for all r. 



E<?o(*HV) > E 50 (*ri V) = AE* 7o (^? c ,r 
>AE* 7 o(*rI*,0 = %o(*riV)- 



) 




(8) 



for all r. 



wherever both quantities are defined. Therefore, by Fubini's theorem, P-a.s., 
we have that (8) holds for £ 
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Provided x is not equidistant from two centers, (8) implies that ^ GS (x) = 
^ GSc (x). Hence, P-a.s., for £-a.e. x, we have ^ GS (x) = ^ GSc (x). Finally, 
using Lemma 15 again, P-a.s., we have, for any stable allocation ip of E = [IT] 
and £-a.e. x, 

g x (^n S ,r)>g x (i;,r)>g x (^ Sc ,r)=g x (^ s ,r) for all r, 

and so we have equality throughout. Thus, any two stable allocations to [II] 
agree for £-a.e. □ 

Finally, in this section we prove Lemmas 15 and 16. The proof of Lemma 16 
proceeds via the following lemma (see [2] and [7] for generalizations). 

Lemma 17 (Mass-transport principle). Let m:'L d x Z d — ► [0, oo] be a 
function satisfying m(u + z,v + z) = m(u, v) for all z € Z d . Then 

m(0, v) = m(u,0). 

v&L d uez d 

Proof. 

m(0,v) = V] m(— v,0) = ^ m(u,0). 

«GZ d v£Z d uGZ d '-' 

For z € Z d , we define the unit cube Q z = [0, l) d + z C R d . 

Proof of Lemma 16. We shall apply the mass-transport principle 
(Lemma 17) to 

m(u,v) = EC{x e Q u -^u{x) € Q v , \x - *n(a;)| < r}. 

Using Fubini's theorem and translation-invariance, we have 

m{0,v)=E£{xeQo:\x-y n {x)\ < r} = P(|0 - ^n(0)| <r), 

v&Z d 

while by a standard property of the Palm process, 

£m(u,0) = E C[^ 1 (OnB(C,r)] = XE*C[^l(0)nB(0,r)]. 

uez d £e[n]nQ 

So the result follows. □ 

The proof of Lemma 15 depends on the two lemmas below. It is convenient 
to prove the first in a general form which will also be useful in proving 
monotonicity results in Section 6. For the present purpose, we may assume 
Ho = S and «o = a - 
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Lemma 18. Let Ho^H be two sets of centers, and let ao > a. Let ijj be 
any stable allocation to 3o with appetite ao, and consider the site-optimal 
Gale-Shapley algorithm for 3 with appetite a. Then for C-a.e. site x, we 
have that tp(x) never rejects x. 

Lemma 19. Let 3 be a set of centers. Let tp be any stable allocation to 
3 and consider the center- optimal Gale-Shapley algorithm for 3. Then for 
C-a.e. site x, we have that x never rejects ip(x). 

Proof of Lemma 18. The proof is by induction on the stage. Suppose 
that the statement is false, and let n be the first stage at which it is violated. 
So at stage n, some center £ G 3o rejects a set of sites R C V ;_1 (£) with 
CR = e > 0. Then at stage n, the center must have shortlisted a set S with 
CS = ao, and with every site in S closer to £ than some site in R is to £. 
Since ip is an allocation, we have £V -1 (0 < a < 

ao, so there must exist a 
set T C S, disjoint from ^~ 1 (£), with CT > e. Since £ shortlisted S at stage 
n, every site x € T must earlier have been rejected by all closer centers in 
3. Hence, by the assumption, for £-a.e. x € T, we have \ip{x) —x\>\£ — x\. 
But this implies that (x,£) is an unstable pair for ip. □ 

Proof of Lemma 19. The proof is by induction on the stage. Suppose 
that the statement is false, and let n be the first stage at which it is violated. 
So at stage n, some nonnull set of sites R C V ;_1 (£) al l reject some center 
£. Therefore, each site in R shortlisted some closer center at stage n, and 
we may therefore find some center £' and some nonnull set T C R such that 
every site in T shortlisted In particular, £' applied to T at stage n. Now 
consider V' _1 (£ / )- ^ every site in is closer to £' than T is, then £' 

must earlier have rejected some nonnull subset of V' _1 (£ / )i which contradicts 
our assumption. Therefore, there exist sites x 6 R' and y € ijj" 1 ^') with 
|y — £'| > |x — £'|. But then is an unstable pair for t/j. □ 

Proof of Lemma 15. The first inequality in (i) is equivalent to the 
assertion that, for C-a.e. x, we have 

(9) \x- V> GS (x)| < \x-i/>(x)\, 

where we write \x — oo\ = oo. Note that in the site-optimal Gale-Shapley 
algorithm, for any x not equidistant from two centers, we have that ip GS (x) 
equals the closest center to x which never rejects x, or oo if all centers reject 
x. Therefore, the first inequality in (i) follows from Lemma 18 with 3o = 3 
and «o = ol. 

Next we consider the second inequality in (ii). Suppose for a contradiction 
that 7^(^,r) < j^(ip GS ,r). Since ^ GS is an allocation, this implies that, in 
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ip, the center £ is unsated or has some site outside B(£,r) in its territory. 
Also, there exists an £-nonnull set T C B(£,r) n [(V' GS )~ 1 (£)\V ,_1 (£)]- From 
(9), for £-a.e. site x G T, we have |cc — V ;GS ( ;z; )| < |s — £|- But now (sc,£) is 
unstable for ip. 

Next we turn to the first inequality in (ii). Suppose, on the contrary, 
that 7^(^ GSc ,r) < 7^(^,r). Since ip GSc is an allocation, this implies that, in 
ip GSc , the center £ is unsated or has some site outside B(£,r) in its territory. 
Also, there exists an £-nonnull set T C B(£,r) n \ (^ GSc ) _1 (0]- But > 

considering the center-optimal Gale-Shapley algorithm, this implies that £ 
rejected T, and this contradicts Lemma 19. 

Finally, we prove the second inequality in (i). Suppose, on the contrary, 
that, for every site x in an £-nonnull set T and two fixed centers we 
have if) GSc (x) = £ and ip(x) = where \x — £| < |x — Then either £ must 
be unsated in ip, or ^ (£) must contain a site y further from £ than a.e. site 
in T, otherwise we would have a contradiction to the first inequality in (ii) 
just proved. But now for a.e. x E T, we have that the pair is unstable 

for □ 

For the remainder of the article, we take \P = \fn, where, Vl/ # is any factor 
allocation, unless otherwise stated. For definiteness, we may take = \P. . 

5. Phases. 

Proof of Theorem 4. Using translation-invariance, the result is an 
immediate corollary of Proposition 20 below. □ 

Define the residual appetite of a center £ to be = = a — 

Proposition 20. Lei II 6e an ergodic point process of intensity A £ 
(0, oo). We /iaue 

P(0 is unclaimed) = (1 — Act) V 0; 

E*£/(0) = (a - A -1 ) V 0. 

Proof. Lemma 16 with r = oo states that P(0 is claimed) = AE*£^' ri J(0). 
Therefore, from the definition of residual appetite, we have 

AE*[/(0) - P(0 is unclaimed) = Xa - 1. 

Recall that no stable allocation may have both unclaimed sites and unsated 
centers. Also, by ergodicity, the existence of unsated centers and the exis- 
tence of an £-nonnull set of unclaimed sites are zero-one events. Therefore, 
at most one of the two terms on the left can be nonzero. The result follows. 

□ 
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6. Monotonicity. We present several results concerning monotonicity in 
stable allocations. Proposition 21 will be used in the proofs of Theorems 6 
and 7. Recall from Section 4 the quantity 

which measures how good an allocation is for a given center. 

The following result states that adding extra centers makes sites happier, 
while making the existing centers less happy. 

Proposition 21 (Monotonicity in H). Fix a, let Hi C H 2 be two sets of 
centers, and let tpi>' l p2 be their respective site-optimal Gale-Shapley alloca- 
tions. 

(i) For C-a.e. site x, we have \x — ip\{x)\ > \x — tp2(x)\. 

(ii) For every center £ £ Hi and all r G [0, oo] , we have 7^ (V>i , r) > 7^ ( V>2 > f) . 

The following result states that increasing the appetite makes sites hap- 
pier, while making centers more unhappy (in an appropriate sense). 

Proposition 22 (Monotonicity in a). Fix H, let a\ < 02, and let ip\,ip2 
be the respective site- optimal Gale-Shapley allocations to H. 

(i) For C-a.e. site x, we have \x — ipi(x)\ > \x — ife(x)\. 

(ii) For every center £ G H and all r £ [0, oo], we have a\ — J^(ipi,r) < 
«2 -7e(^2,r). 

Proof of Proposition 21. (i) We apply Lemma 18 with H = Hi, 
H = H2, tp = ipi, ip ='4 ) 2, a$ = a. The result follows upon recalling that 
the site-optimal Gale-Shapley algorithm allocates a site to the closest center 
which does not reject it. 

(ii) Suppose the statement fails for some £ and r. Then there exists a 
positive volume of sites x € B(£,r) with ^(x) = £ but ipi(x) {£, A}, and, 
furthermore, £ covets all such sites under ipi. Pick such an x which is not 
equidistant from any two centers. By Proposition 21 (i), we have \x — i/ji(x)\ > 
\x — £|, and since x is not equidistant from any two centers, the last inequality 
is strict. But now x desires £ under ipi, so (x,£) is an unstable pair for ijji. 
□ 

Proof of Proposition 22. (i) We apply Lemma 18 with H = H, 
a = a\, ao = OLii tp = ipi, ip GS = ip2- The result follows on recalling that the 
site-optimal Gale-Shapley algorithm allocates a site to the closest center 
which does not reject it. 

(ii) Suppose the statement fails for some £ and r. Then there exists a 
positive volume of sites x G -B(£, r) with ^(x) = £ but ipi(x) {£, A}, and, 
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Fig. 3. An illustration of the proof of Theorem 5(i). If the all centers illustrated have 
unbounded territory, we obtain a contradiction to stability. Here form an unstable 

pair. 

furthermore, £ covets all such sites under ipi. Pick such an x which is not 
equidistant from any two centers. By Proposition 22 (i), we have \x — \j}\(x)\ > 
\x — £|, and since x is not equidistant from any two centers, the last inequality 
is strict. But now x desires £ under ipi, so is an unstable pair for ip±. 

□ 

7. Geometry of territories. In this section we prove Theorem 5. State- 
ment (iii) of that theorem will follow from a somewhat stronger result, The- 
orem 24, in which we construct a "canonical" version of the stable allocation 
which has several good properties. 

A key tool in proving Theorems 5 and 6 is the following geometric lemma 
inspired by arguments in [3]. Figure 3 illustrates the main idea. Let 5 = {i£ 
IR 6 ': \x\ = 1} be the unit sphere. A cap is a proper subset of S of the form 
H = S n B(y, r), where y € S. A cone is a set of the form 

V = Vh = {ah :h£H and a 6 (0, oo)}, 

where if is a cap. 

Lemma 23 (cones). For each d>l, there exist an integer k and cones 
V±, . . . , Vfe whose union is M. d \ {0}, such that, for any x, y lying in the same 
cone Vi, 

if \x\ > \y\, then \x — y\ < \x\ — \\y\. 
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(For example, when d = 2, we may take k = 13 and Vi = {x ^ : argx E 
[27r(i - l)/13,27r</13)}.) 

Proof. The set of all caps of diameter 1/2 is a cover of the compact 
set S, hence, it has a finite sub-cover H\,... say. Put Vi = Vh v 

Suppose x,y E Vi and |x| > |y|. Let f3 = \y\/\x\, so that \(3x\ = \y\, and 
hence, f3x,y E Then we have 

\x — y\ < \x — (3x\ + \(3x — y\ 

< (1 - + diam(\y\Hi) = \x\ - \y\ +\\y\. □ 

Proof of Theorem 5. (i) Without loss of generality, we may assume 
that IT is ergodic under translations; if not, we may appeal to the ergodic 
decomposition theorem ([11], Theorem 10.26). 

Call a center bad if its territory is unbounded. Supposing for a contradic- 
tion that bad centers exist, they must form an ergodic translation-invariant 
point process of positive intensity. Let V±, . . . , Vk be cones as in Lemma 23, 
and note that, since any cone contains arbitrarily large balls, it contains 
infinitely many bad centers a.s. For R > 10, define the sets 

To = 5(0,1); 

Ti = Vi n [B(0, R) \ B{0, 10)], i = l,...,k. 

The above observations imply that if R is chosen sufficiently large, with 
positive probability, we have 

(10) each of To, . . . , contains a bad center. 

Suppose the event in (10) occurs, and let be a bad center in Ti, for 
each i. Since £o is bad, there exists a site x € M. d \ B(0,R) in its territory. 
And x lies in some cone, Vj, say. Since £j E Vj also, and \x\ > R > > 10, 
Lemma 23 applies to give 

\x — £j| < |x| — g 

< \x\ — 5 

< |a;-^o|- 

Furthermore, since is bad, there exists a site y in its territory with \y — 
£j \ > \x — But then the pair (x,£j) is unstable, a contradiction. 

(ii) Without loss of generality, LT is ergodic. By translation-invariance, 
it suffices to prove that B(0, 1) intersects only finitely many territories a.s. 
Note also that, since H is almost surely discrete, if a ball B intersects in- 
finitely many territories, then ^f(B) must contain centers at arbitrarily large 
distances from B. 
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Consider the integer lattice 7L d . Call a lattice site z € Z d bad if B(z, 1) 
intersects infinitely many territories, and suppose for a contradiction that 
P(0 is bad) > 0. Defining To,...,T^ as in the proof of part (i), a similar 
argument shows that R may be chosen such that, with positive probability, 

(11) each of To, . . . , Tfc contains a bad lattice site. 

On the event in (11), there exists a site x € B(0,2) with € R d \B(R), 
and ^(x) € Vj, say. Let zj be a bad lattice site in Tj\ then there exists 
y € 1) with \y — ^(y)\ > \y — ^(x)\. Lemma 23 gives 

\y-V(x)\ < \zj - + 1 

<|*(x)|-i|z,-| + l 

< -4 

< |* (a;) - x\. 
So the pair (y,^(x)) is unstable. □ 

The statement in Theorem 5(iii) does not hold for an arbitrary stable 
factor allocation, since, for instance, we allow ^(x) = A on an arbitrary C- 
null set. Indeed, it does not necessarily hold for \]/ GS , since $ GS (x) = A for 
all sites which are equidistant from two or more centers. We shall define a 
"canonical version" ^ of the stable allocation for which Theorem 5(iii) and 
other desirable properties hold. 

Let V be a stable allocation, and define a new function ip as follows. 
For Borel sets i,B£ M d , we say A is essentially a subset of B, and write 
A Q css B, if C{B \ A) = 0. For a site x, if there exists £ € H U {oo} such that 
x has some neighborhood which is essentially a subset of tp -1 ^), then we 
put tp(x) = £. If there is no such £, we put ip{x) = A. Note that ip is well 
defined because no neighborhood can be essentially a subset of two disjoint 
sets. Note also that any two allocations which agree £-a.e. yield the same 

PROOF of Theorem 5(iii). Immediate from Theorem 24 below. □ 

For the remainder of the section, we assume that LT is ergodic. 

Theorem 24 (Canonical allocation). Let ^ 9 be an isometry-equivariant 
factor allocation to IT (e.g., ^ GS j and define ^. by the above procedure. Then 
P-a.s., is the unique minimizer of the set ^ _1 (A) in the class of stable 
allocations ip of [IT] having all territories open and the unclaimed set open. 
Furthermore, ^, is an isometry-equivariant factor allocation, and a.s. each 
territory of ^ has only finitely many connected components, and has C-null 
boundary. 
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In order to prove Theorem 24, we first need to establish some properties 
of For a center £, define 

%)=sup{|x-e|:xG^ 1 (0}. 

[Note that may differ between allocations which agree a.e.] Define the 
random set of sites 

Z = \J{x€R d :\x-£\ =R(C)} 

U \J{x€R d :\x-£\ = \x-£% 

where the two unions are over all centers and all pairs of centers respectively. 
Since H countable, Z is £-nuil a.s. 

Lemma 25. P-a.s., for every site x ^ Z , x has a neighborhood N such 
that #(*(JV)\A) = 1. 

Before proving Lemma 25, we make some further definitions. Let B be 
a ball. By Theorem 5(ii), a.s., B intersects only finitely many territories, 
belonging to centers £i, .. say. Define the set 

Z B = \J{x€R d :\x 

i 

u\J{xeR d :\x-&\ = \x-£ j \}. 

Define a cell (of B) to be any connected component of B \ Zb- Since B \ Zb 
is open, all cells are open sets, and the following lemma implies that, for any 
given B, there are only finitely many cells. 

Lemma 26 (Intersecting hyperplanes). Let Y be a union of finitely many 
(d — l)-spheres and (d— 1)- dimensional hyperplanes in R d . Then R d \Y has 
only finitely many connected components. 

Proof. By applying an inversion about a site not in Y, we may assume 
that Y is a union of spheres only, and is therefore bounded. Supposing that 
has infinitely many components, a repeated bisection argument shows 
that there must exist a site z all of whose neighborhoods intersect infinitely 
many components. But considering a neighborhood of z small enough to 
avoid all those spheres which do not pass through z, we see that this is 
absurd. □ 

Lemma 27. P-a.s., if C is any cell, then |^(C) \ A| = 1. 



26 



C. HOFFMAN, A. E. HOLROYD AND Y. PERES 



Proof. We claim first that, for each Zb with ^(x) 7^ A, we have 

that, ^(x) equals the unique closest center to x in the set 

r(z) = {&:|x 

or \J/(x) = 00 if T{x) is empty. To prove this, note first that, by the definitions 
of and Zb, we must have ^(x) £ T(x) U {00}. But if T(x) contains a 
center £j which is closer to x than \I/(x) is, then the pair (x,£j) is unstable. 
This proves the claim. 

Now note that the k 2 quantities 

[a; -&| -.R(&), i = l,...,k, 

I x ~ Ci I — I x ~ Cj I j * 7^ J) 

are continuous in x, and nonzero except on Zb. It follows that their signs 
are all constant on any given cell. The result now follows from the above 
claim. □ 

PROOF of Lemma 25. Let B = B(x,l), and define Zb and cells as 
above. Clearly, we have Zb Q Z, so x lies in some cell. Since all cells are 
open, x has a neighborhood lying in the cell. Now we use Lemma 27. □ 

Proof of Theorem 24. First note that, since Z is £-null, Lemma 25 
implies that a.s. VP agrees with ^ £-a.e. It follows immediately that is an 
allocation a.s. By the construction, inherits the isometry-equivariance of 
Vl/. . It also follows immediately from the definition of ^ that the territories 
and the unclaimed set are open. 

To prove that ^ is a stable allocation, suppose, on the contrary, that 
(x,£) is an unstable pair. Suppose first that 'I'(x) = ( and ^(y) = £, where 
\x — £| < |x — £| and \x — £| < \y — £|. Then x and y have neighborhoods 
intersecting ^~ 1 (£) and l I / ~ 1 (£) in full measure respectively, so we may find 
x',y' with \I/(a/) = ( and = £, and close enough to x,y that we still 

have \x' — £| < \x' — £| and \x' — £| < |y' — £|. Similarly, if x is unclaimed in 

then we may find a nearby x' which is unclaimed in and finally, if £ 
is unsated in ^, then it is unsated in 'I'. In each case we deduce that (x',£) 
is an unstable pair for ^, a contradiction. 

To prove the required minimality statement, we note that a.s., for any tp 
as described in the statement of the theorem and any x with ip(x) 7^ A, we 
have that tp takes the value ip{x) on some neighborhood of x, so ^f(x) = ip{x). 

To prove that the territories of ^ have only finitely many components, 
note that, in view of Theorem 5(i) and translation- invariance, it is sufficient 
to rule out the possibility that the intersection of -B(0, 1) with some territory 
has infinitely many components. Let B = B(0, 1) and define Zb and cells 
as above. Recall that there are only finitely many cells (by Lemma 26), 
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and that each cell is connected. Lemma 27 combined with the definition 
implies that is constant on each cell. Hence, it is sufficient to rule out the 
possibility that the intersection of some territory with B has a component 
lying entirely in B n Zb- But this is impossible since any such component is 
open, while Zb is £-null. 

Finally, to prove that the territories of ^ have £-null boundaries, note 
that, by Lemma 25, if x ^ Z, then x does not lie in the boundary of any 
territory. But Z is £-null a.s. □ 

8. Infinite demand. In this section we prove Theorem 6. We shall deduce 
it from a stronger result which applies to more general point processes. 

For a simple point process LT and a Borel set S C M. d with CS € (0, oo), 
let II s be the point process obtained by adding a random extra center in S: 

U s (A)=U(A)+5u(A), 

where U is uniformly distributed in S, independent of LT. We call the process 
LT insertion-tolerant if the law of H s is absolutely continuous with respect to 
the law of LT for every S. Let LT5 be the point process obtained by deleting 
all centers in S: 

U S (A) = U(A\S). 

We call LT deletion-tolerant if the law of Us is absolutely continuous with 
respect to the law of LT for every 5. 

Theorem 28. Let U be an ergodic point process and let \a = 1. 

(i) //LT is insertion-tolerant, then a.s. every center is desired by an in- 
finite volume of sites. 

(ii) If IT is deletion-tolerant, then a.s. every site is coveted by infinitely 
many centers. 

Proof of Theorem 6. It is easy to check that a Poisson process is both 
insertion-tolerant and deletion-tolerant, so the result follows immediately 
from Theorem 28. □ 

Next we prove Theorem 28. Our starting point is the following weaker re- 
sult which states that an infinite set of sites would desire some appropriately- 
positioned center in the unit ball. This result will also be used in the proof 
of Theorem 7. In the following, we take \x — 00 1 = 00 for x £ W 1 . 

Lemma 29. Let Xa = 1 and let LT be an ergodic insertion-tolerant pro- 
cess. We have a.s. 

£{xeR d :\x-^(x)\ > \x\ -l} = oo. 
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V - 2u V - u 



Fig. 4. An illustration of the proof of Theorem, 28(i). Some cone must contain an infinite 
volume of sites which desire all centers in the ball at its apex, and therefore, it must also 
contain an infinite volume of sites which desire all centers in the ball B(0, 1). 

Proof. Fix K < oo. We shall show that the quantity in question is at 
least K. Let n' be a point process obtained from LT by superposing m extra 
centers, independently uniformly distributed in the ball B = B(0, 1), where 
m is an integer such that am > K. By repeatedly applying the insertion- 
tolerance condition, the law P' of LT' is absolutely continuous with respect 
to the law P of 11. Hence, if is a stable factor allocation as usual, Vl/n' is 
defined C-a.e. P'-a.s. Also, since Xa = 1, by Theorem 4, P-a.s., every center 
is sated in v&n- Hence, the same holds P'-a.s. in ^n'- 

The last observation implies that, P'-a.s., we have 



Now since [II] C [n'], Proposition 21 (i) gives that, almost surely with respect 
to the joint law, for £-a.e. x € ^^(B), we have 

[a; - *n(a;)| > \x - ^n'(^)| > M - 1- 
Combining this with (12) yields that, P-a.s., 



Proof of Theorem 28. (i) The proof is illustrated in Figure 4. Let 
V±, . . . , Vk be the cones of Lemma 23. Since they partition M rf \{0}, Lemma 29 
implies that, for some fixed £, we have 



(12) 



C^ u }(B)>am>K. 



C{x G R d : \x - Vn(x) \ > \x\ -1}>K. 



□ 



(13) 



P[£{x£V e :\x-^(x)\ > \x 



l} = oo] >0. 
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Pick such an t and write V = Vg, and also pick w € M rf such that B(w,l) C 
V \ B(0, 10). Note that, for any n > 1, we also have 

(14) B(mi;,l) CV\S(0,10). 

For z £ M d and ^4 C we denote the translated set A + z = {a + z:a£ 
Define the set of sites 

L n = {x £ V — nw : \x — ^(x)\ > \x + nw\ — 1}, 

and the event 

A n = {C{L n ) = oo}. 

The sequence of events (i n )n>o is stationary and ergodic under P, and (13) 
implies that P(^4q) > 0. Therefore, we have 



(15) 



1. 



.71=1 

Now suppose that the event A m occurs, where m > 1. We shall deduce 
that every center in B(0, 1) is desired to be an infinite volume of sites. Let 

z G L m \ B(—mw, m\w\ + 1) 

(the occurrence of A m guarantees that an infinite volume of such sites z ex- 
ists). Let £ be any center in B(0, 1) (if one exists). We shall apply Lemma 23 
with x = z + mw and y = £ + mra. Note that z £ V — tow implies z + tow S V, 
while (14) implies £ + mw £ V. Furthermore, since z £ B(—mw,m\w\ + 1), 
we have \z + mw\ > m\w\ + 1 > \mw\ + |£| > |£ + mw\. Thus, Lemma 23 
applies to yield 

\z — £| < \z + mu>| — + mw|. 

Since z 6 L m , we have |z + tow| < \z — \&(z)| + 1. By (14), we have |£ + m«;| > 
10. Substituting into the above equation, we obtain 

\z-£\<\z- *(z)| + 1 - 5 < \z - 

Hence, 2 desires £. 

The above argument together with (15) shows that, almost surely, every 
center in B(0, 1) is desired by an infinite volume of sites. Since M. d may be 
covered by a countable collection of unit balls, by translation invariance, the 
same claim applies to every center. □ 

Now we turn to the proof of Theorem 28 (ii), which follows similar lines 
to that of (i). For a center £, define the radius of its territory: 

R((,) = = esssup |£ — x\. 
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Lemma 30. Let Xa = 1 and let U be an ergodic deletion-tolerant process. 
We have a.s. 

#{£GE :j R(£)>|£|-l} = oo. 

Proof. We mimic the proof of Lemma 29. We show that the quantity 
in question is at least K. Write B = B(0, 1). By scaling M. d , without loss of 
generality, we may assume that aK < CB, and the process is still deletion- 
tolerant. Let n / be the process obtained by deleting all centers in B. By 
Theorem 4 and deletion-tolerance, for II' a.s., infinitely many centers desire 
some site in B. Hence, by Proposition 21(h), the same is true for II. □ 

Proof of Theorem 28. (ii) We mimic the proof of (i). Take a cone V 
such that 

P(#{£ G V n S : R(0 > 1^1 - 1} = oo) > 0, 

and take w as before. Almost surely some event 

A n = {#{£ G (V - nw) n S : > |f + nw\ - 1} = oo} 

occurs, and the previous argument goes through to show that infinitely many 
centers desire every site in B(0, 1). Therefore, the same is true for every site. 
□ 

9. Critical lower bounds. In this section we deduce Theorem 7 from more 
general results. The following is Corollary 15 in [10], where it was proved 
making use of results from [9] and [13]. 

Theorem 31. For U a Poisson process with A = a = 1 and for ^, any 
(not necessarily stable) factor allocation with a.e. site claimed and every 
center sated, we have EX d / 2 = oo for d= 1,2. 

Interestingly, for d > 3, there exist "transport rules" (these are a mild gen- 
eralization of factor allocations) with exponential tails; see [10] for details. 
We shall prove the following. 

Theorem 32. Let II be ergodic and insertion-tolerant and let Xa = 1. 
For any d>\, we have EX rf = oo. 

Proof of Theorem 7. Immediate from Theorems 31 and 32 above 
[and Theorem 4(h)]. □ 

We prove a somewhat weaker statement for deletion-tolerant processes. 
Recall the definition of the radius R(£) from the last section, and write 
R*(C) = Ry u * (£) for the radius associated with the Palm process. 
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Theorem 33. Let II be ergodic and deletion-tolerant and let Xa = 1. 
For any d>l, for the Palm process II*, we have E*R*(0) d = oo. 

Proof of Theorem 32. Since II is insertion-tolerant, Lemma 29 ap- 
plies. Taking the expectation of the quantity in Lemma 29 and using Fubini's 
theorem and translation-invariance, we have 

BC{x:\x - $>(x)\ > \x\ - 1} 

= / P(|x-*(x)|>|x|-l)£(dx) 

= f P(X > \x\ - l)C(dx). 
For a constant c = c(d) £ (0, co), the last integral may be written 

pod 

/ P(X >r- l)cr d_1 dr. 

Jr=0 

Writing \x for the law of the random variable X + 1 and using Fubini's 
theorem again, the above equals 

c / r d ~ l dr fi(dz) = -E[{X + l) d }. 
Jz=0Jr=0 d 

Hence, Lemma 29 implies that the above quantity is oo. □ 

Proof of Theorem 33. The expectation of the quantity in Lemma 30 
equals 

/ P*(R*{0)>\x\-l)£(dx); 

JM. d 

so this quantity is infinite. We now proceed as in the proof of Theorem 32. 

□ 

Open problems. 

(i) Nonuniqueness. Does there exist a discrete set of centers in M. d which 
admits two stable allocations (differing on a nonnull set)? By Theorems 2 
and 3, such a set must have greater than quadratic growth, and may not 
arise from a translation- invariant point process. 

(ii) Critical tail behavior. What is the tail behavior of X for the critical 
Poisson model? In particular, give any quantitative upper bound for d>2. 

(iii) Percolation. In the Poisson model, does the set of unclaimed sites 
have an infinite component for a sufficiently small? Does it fail to have an 
infinite component for sufficiently large a < 1? Similar questions apply to the 
set of claimed sites. After this paper was accepted for publication, progress 
on these problems was achieved in [5]. 
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(iv) Connected territories. Is there a translation-equivariant allocation 
(not stable) of the Poisson process in which every territory is connected, in 
the critical two-dimensional case (d = 2 and Xa = 1)? 

(v) Hyperbolic space. The concept of a stable allocation extends to more 
general settings, such as hyperbolic space in place of M. d . The proofs of 
Theorems 1, 3 and 4 adapt to this setting, but the proof of Theorem 5 does 
not. Are all territories bounded in the setting of hyperbolic space? 
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